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Abstract 

Motivated by the design of an integrated CMOS-based detection platform, a simulation model for CCD and CMOS imager-based lumines- 
cence detection systems is developed. The model comprises four parts. The first portion models the process of photon flux generation from 
luminescence probes using ATP-based and luciferase label-based assay kinetics. An optics simulator is then used to compute the incident 
photon flux on the imaging plane for a given photon flux and system geometry. Subsequently, the output image is computed using a detailed 
imaging sensor model that accounts for photodetector spectral response, dark current, conversion gain, and various noise sources. Finally, 
signal processing algorithms are applied to the image to enhance detection reliability and hence increase the overall system throughput. To 
validate the model, simulation results are compared to experimental results obtained from a CCD-based system that was built to emulate the 
integrated CMOS-based platform. 
© 2004 Elsevier B.V. All rights reserved. 
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1. Introduction 

Conventional biological assays are highly repetitive, la- 
bor intensive, and require microliter volume samples. The 
associated biochemical protocols often require days or 
weeks to perform at a cost of hundreds of dollars per test. 
Problems remain in detecting and quantifying low levels 
of biological compounds reliably, conveniently, safely and 
quickly. There is also a growing interest in the development 
of inexpensive portable biosensors for environmental and 
biomedical diagnostics. Solving these problems will require 
the development of new techniques and sensors, not only 
to selectively identify target compounds, but also to assay 
large numbers of samples. 

A biochemical testing procedure can be divided into four 
steps: sample preparation, assay, detection, and data analysis 
as shown in Fig. 1 . Currently, each step is being separately 
automated and miniaturized. However, there continues to 
be a need for designs that accommodate efficient integrated 
circuit manufacturing techniques to realize associated cost 
savings. We have been investigating the integration of three 
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of these four steps into a single miniaturized platform as 
shown Fig. 2. 

A variety of assay methods have been developed for 
molecular detection. These methods include electrochem- 
istry (Woolley et al., 1998), optical absorption (Kunz, 1997), 
interferometry (Verpoorte et al, 1992), luminescence and 
fluorescence (Haugland, 1998). In this paper we focused on 
luminescence detection or luminometry. This technique is 
becoming increasingly popular due to its high sensitivity, 
wide dynamic range, and relatively inexpensive instrumen- 
tation. Superior sensitivity and low background distinguish 
luminometry from other analytical methods. Luminometry 
is up to five orders of magnitude more sensitive than absorp- 
tion spectroscopy and more than 1000 times more sensitive 
than fluorometry. A state-of-the-art luminometer can detect 
as little as 0.6 pg of adenosine triphosphate (ATP) or 0.1 fg 
of luciferase (M 100 molecules), two common luminescent 
analytes (Turner et al, 1985). Numerous bioluminescent 
and chemiluminescent reactions are studied using lumi- 
nometry and are commonly used in biotechnology research, 
environmental testing, industrial applications, and clinical 
research. Among its many applications are the measure- 
ments of gene expression using reporter gene assays, the 
determination of intracellular ATP, and DNA sequencing. 

Commercially available platforms for photon detection 
use bulky CCD camera-based setups and require the use 
of large quantities of reagents due to the light loss in the 
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Fig. 2. Proposed system for future biochemical testing. 

optical path. As a result these systems are not well suited 
to low cost applications. The development of imaging sen- 
sors in standard complementary metal-oxide-semiconductor 
(CMOS) technologies makes it possible to integrate sensing 
and processing on the same integrated circuit, enabling many 
low power and low cost applications. We have designed an 
integrated CMOS-based detection platform for use with lu- 
minescent microarrays. (Eltoukhy et al, 2004) A major con- 
cern in the design of such a miniaturized system is that the 
photogenerated signals can be very weak and therefore dif- 
ficult to detect reliably. To quantify the detection limits of 
the envisioned system for the required assay concentration 
and throughput, we built a complete model for simulating 
the path from photogeneration through the optical path to 
detection and image post processing. Such modeling is used 
to guide the overall system design, whereby the special char- 
acteristics of a wide array of luminescent assaying methods 
can be exploited to enhance detection sensitivity beyond that 
of off-the-shelf CCD or CMOS-based sensors. 

In Section 2, we describe the light generation process 
of luminescent probes and formulate the complete kinetic 
model of ATP-based and luciferase label-based assays. In 
Section 3, we describe the model of the proposed lumi- 
nescent detection system including the associated optical 
pathway, imaging array characteristics, and applied signal 
processing algorithms. Finally, in Section 4, we present 
simulation and experimental results. 



2. Luminescence light generation processes 

Luminescence assay techniques are divided into two gen- 
eral categories. The first is direct target detection in which 
the photon emitting species physically interacts with the 
target of interest at a predetermined location. An example of 
this approach is luminescence-based immunoassays, which 
usually involve probing a protein of interest with a primary 
antibody that reacts with a secondary antibody. Light is pro- 
duced when the enzyme bound to the secondary antibody 



acts upon a luminescent probe. The second category is in- 
direct or linked detection in which the luminescent species 
indirectly measures the targeted characteristic, usually 
through an intermediate chemical process (Tang et al., 1995; 
Koster et al., 1997; Van Dyke et al, 2002; De Mello, 1996). 
While multiplexing is feasible in direct target detection, in- 
direct detection requires confinement of the photon genera- 
tion process as well as physical barriers for the independent 
reactions. This physical isolation complicates the integra- 
tion of indirect detection methods to miniaturized systems, 
since high density mechanical barriers (e.g. micro-wells) 
as well as micro-scale solution delivery systems would be 
required. Since indirect detection is the more general of the 
two categories, it is the focus of the following analysis. 

2.1. Luminescence light generation 

The time-dependent light generation from a typical lumi- 
nescence process is a function of the underlying chemical 
reaction kinetics. The rate of a reaction, in general, is the 
speed at which reactants are converted into products. If en- 
zyme species E (the catalyst) converts the substrate molecule 
S into product P 9 the stoichiometric formula is given by 



S+E kt ^ T E+ P, 



(1) 



where kf and k T are the association and disassociation rate 
constants. In (1) the reaction rate is denned by 



rate = 



d[P] d[S] 



dt 



dt 



= kf[S][E]-k T [E][Pl 



(2) 



where [E], [S], and [P] are the concentrations of the enzyme, 
substrate and product in the medium, respectively. Now if we 
suppose that the above process is a luminescence enzymatic 
reaction with quantum yield a, then the photon generation 
rate / in volume Vof the reaction medium (A is Avogadro's 
number) would 'be 

/ = («^)~~ = (aVA)[E](k f [S) - k r [P]). (3) 

The total number of photons generated by this lumines- 
cence process iVph (7), in the time interval T, would be 



(4) 



N ph (T) = (aVA)[E]j (k f [S] - k T [P])dt 

In luminescence assays the experiment is typically set up 
in such a way that the luminescence probe (e.g. a light gen- 
erating enzyme) either reports the quantity of a substrate 
molecule (e.g. ATP) or the molecule to which it binds (e.g. 
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luciferase-based labels in immunoassays). The photon gen- 
eration rate from the luminescence reaction, which is a func- 
tion of the target concentration, is then measured and corre- 
lated to the target concentration. 

2.1.1. Substrate detection kinetics 

In the first group of luminescence assays, the rate at which 
photons are generated represents the substrate concentration 
given by (3). As the substrate is consumed by the catalyst, 
the light intensity decreases and eventually approaches zero. 
If we assume that the disassociation rate is insignificant (i.e. 
negligible inhibition), the light intensity with initial substrate 
concentration [So] becomes 



7(f) = {aVA)^^ = (aVA)kf[E][S 0 ]e- kf[E]t . 
at 



(5) 



The total amount of photons from time t = 0, the start of 
the process, to t = T is 

N ph (T) = (aVA)[S 0 ]{\ - e"*^ r ). (6) 

The photon intensity in (6) is proportional to the target 
concentration, but the light intensity decays exponentially 
with a time constant, which is a function of the catalyst 
concentration and turnover rate k{ . 

2.1.2. Catalyst detection kinetics 

The second approach in luminescence assays is to link 
the target molecule quantity to a luminescence catalyst. In 
this approach, excess substrate is used, making sure that 
its consumption does not affect the reaction kinetics. If the 
saturation concentration for the substrate is [S max ] and the 
target concentration is equal to the catalyst [E\, the light 
intensity becomes 

/(,) = {aV A)^p- = (aVA)k f [S max ][E]. (7) 

Since the light intensity based on (7) is time independent 
and proportional to the target (or the catalyst) concentration, 
the total number of photons generated from this process is 
a function of the integration time, i.e. 



N ph = (aVA)k { [Sm^][E]T. 



(8) 



In such assays, the number of target molecules can be very 
small resulting in a low but steady light. Long integration 



times, however, can be used to collect a significant number 
of photons. 

2.2. Models of specific luminescent assays 

2.2.1. ATP measurement 

ATP bioluminescence assays are designed to measure the 
quantity of adenosine 5 '-triphosphate (ATP) in a sample 
(Schena et al, 1995). Its applications include indirect mea- 
surement of bacteria, yeasts, fungi and other microorgan- 
isms, which have a regulated number of ATP in foodstuffs, 
beverages, water and other media. The assay typically 
employs recombinant luciferase to catalyze the reaction 



ATP + luciferin + O2 luciferas > e AMP + oxyluciferin 

+C0 2 + PPi + /zv. 



(9) 



In most practical assays, the concentration of luciferin is 
high enough that we can consider it to be in deep saturation 
as shown in Fig. 3a. In this case, and when product inhibition 
is negligible, the rate of ATP consumption is given by 

= -£ L [£][luciferin^ 

(10) 

where k^ is the association rate of luciferase macro-reaction. 
We can rewrite (10) by substituting k t — fc L [£][luciferin max ] 
as the turnover rate of the ATP consumption to obtain 



7(r) = (aVA)k t [AT? 0 ]Q- ktt . 



(11) 



The quantum efficiency of luciferase is about 0.88, and 
the turnover rate of the enzyme, depending on the ratio of 
luciferase to ATP, can vary from 0.1 to Is -1 (glow com- 
pared to flash). In the glow process, the light intensity is 
approximated by 

7(f) « 5.3 x 10 22 V[ATP 0 ]e~ ou (12) 

2.2.2. Luminescent labels 

If the reporter of a specific biological process is the cata- 
lyst of a luminescence reaction, the assay can be optimized 
in such a way that the rate-limiting factor becomes the cat- 
alyst concentration (Kricka, 1988). In such an assay, any 
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(a) (b) 
Fig. 3. Block diagram of (a) ATP detection assay and (b) luciferase detection system. The (00) symbol indicates excess concentration. 
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change in target concentration changes the catalyst concen- 
tration, therefore altering the photon flux intensity as shown 
in Fig. 3b. The light intensity in the case where luciferase 
is the label of the target becomes 

/(f)* 5.3 x 10 22 V[£ 0 ]. (13) 

223. Pyrosequencing 

Pyrosequencing is a DNA sequencing method based on 
the detection of released inorganic pyrophosphate during 
DNA synthesis. Using a linked enzymatic reaction, visible 
light (at 562 nm) proportional to the number of incorpo- 
rated nucleotides is generated (Brovko et al., 1994; Ronaghi 
et al., 1996; Ronaghi, 2001). The enzymatic cascade begins 
with a DNA polymerization reaction in which inorganic py- 
rophosphate (PPi) is released as a result of nucleotide incor- 
poration by polymerase. The released PPi is subsequently 
converted to ATP by ATP-sulfurylase. The synthesized 
ATP provides the energy for luciferase to generate photons. 
Unincorporated deoxy-nucleotides and ATP are degraded 
by the enzyme apyrase to chemically reset the enzymatic 
system after the incorporation test. The enzymatic reactions 
in this method are given by 

(DNA), + dNTP ?0lymeraS > e (DNA) n +i + PPi (14) 



A ^ ATP-Sulfury lase irT _ _ 9 

PPi + APS ATP + SO4- 2 



ATP -I- luciferin + O2 luciferas > e AMP + oxyluciferin 



(15) 



dNTP dNMP + 2Pi, 



+C0 2 + PPi + /zv (16) 
(17) 



where APS is adenosine phosphosulfate, AMP, adenosine 
monophosphate, dNTP, deoxynucleotide triphosphate, and 
Pi, phosphate. As shown in Fig. 4 this enzymatic system 
regulates the light generation by recycling PPi, but the de- 
grading enzyme, apyrase, breaks all nucleotide and ATP 
molecules in time, thus decreasing the light intensity. If we 
assume that PPi regulation has a unity gain positive feed- 
back, and the turnover rate of apyrase is & a , then the light 



generated by single incorporation from this bioluminomet- 
ric assay is 

7(0 = (aVA)k t [DNA]Q~ kat . (18) 

For most practical applications, & a ^ 0.05 s" 1 and k t & 
Is -1 , and thus the light intensity becomes 

I{t) = 5.3 x 10 23 V[DNA]e-°' 05/ . (19) 

In (19) the negative effects of product inhibition are not 
included. While these unwanted pathways can potentially 
alter the kinetics of the reaction and the light intensity in 
general, one can still approximate light intensity in cases in 
which the sample concentration is low. 



3. Detection system modeling 

Having established a comprehensive model for the quan- 
tum efficiency of bioluminometric assays, we now discuss 
the detection portion of the system. Various detection 
systems have been developed for luminometry. By far 
the most sensitive detection devices are photomultiplier 
tubes (PMTs), which via a photocathode and a series of 
amplifying dynodes can generate up to one million elec- 
trons for every incident photon. These devices are very 
sensitive since their noise can be removed easily using a 
level-discriminator. However, they are costly and require 
high operating voltages (1000 VDC), precluding their use 
in a portable system. Furthermore, the overall photon detec- 
tion efficiency of PMT-based systems is limited to 1-4% by 
the optics and the low quantum efficiency (10%) of PMTs. 
Finally, multiplexed imaging is not readily feasible using 
PMT's, as they are relatively large and thus unsuitable for 
dense arrays. It is this aspect of PMT's that has most limited 
their use and applicability (Van Dyke et al., 2002). 

The most commonly used visible range imaging sensor 
architectures employ either charge-coupled device (CCD) 
or CMOS photosensor arrays. The major difference be- 
tween the two is the specific readout mechanism used. 
CCDs employ a "bucket brigade" to serially shift out the 
photogenerated electrons accumulated at each photosite. 



dNTP 



Target Molecule 

[DNA] \ [A PS] 

J — . V L 



[Polymerase] 




[Luciferin] 



\i7 



[Apyrase] 



[ATP- 
Sulfurylase] 



ATP 



[Luciferase] 



PPI h\ 



[Apyrase] 



Fig. 4. Block diagram of pyrosequencing in which PPi released from DNA polymerization is measured. The (00) symbol corresponds to excess 
concentration. Apyrase competes with polymerase and luciferase for dNTP and ATP consumption. 
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CCDs are fabricated in a nonstandard semiconductor pro- 
cess that is solely optimized for sensing and charge transfer 
(Hoist, 1991). As a result, CCD image sensors achieve very 
high sensitivity, low noise, and low non-uniformity; hence 
they have earned a central place in the biological imaging 
arena. Although PMTs are more sensitive, unfortunately it 
has been our experience that ensuring the reliability of the 
chemistry itself places the true lower bound on detection re- 
quirements and thus CCDs do not have practical deficiencies 
in this respect. Since the thermal and shot noise generated 
from the photodiode junction cannot be differentiated from 
the photogenerated signal, cooling the device is the primary 
means of reducing such noise to negligible levels. Although 
liquid nitrogen is used in applications requiring extreme 
sensitivity (77 K), a lower cost and typically adequate so- 
lution involves incorporating a stack of Peltiers to cool the 
sensor down to as low as 200 K. This, however, requires the 
use of several high supply voltages resulting in high power 
consumption. Moreover, no other analog or digital circuits, 
such as for clock generation, timing, analog-to-digital (A/D) 
conversion, digital processing and storage, can be inte- 
grated with a CCD image sensor on a single chip resulting 
in multi-chip imaging system implementations with high 
power consumption, high cost, and large size. 

CMOS, on the other hand, offers the benefits of CCDs 
with the possibilities of circumventing many of their draw- 
backs. For instance, CMOS offers complete customizability 
of the photodetection array to suit the requirements of the 
specific biological application. Hence the detection area 
of each photodiode can be optimally sized with respect to 
the assay volume and its light generation characteristics. 
Furthermore, this quality allows one to forego intermedi- 
ary optics and to perform "contact" photonic detection of 
chemiluminescent assay signals, thereby drastically increas- 
ing photodetection efficiency to near QE limited percentages 
while eliminating expensive and bulky optics. This improve- 
ment alone can increase sensitivity of the system by an order 
of magnitude. Although others (Eggers et al., 1994; Lamture 
et al., 1994) performed such detection using a CCD-based 
sensor, it is the leveraged use of an integrated CMOS process 



that leads to considerable overall gains in system perfor- 
mance. More specifically, CMOS makes true low power op- 
eration possible, with the ability to integrate both the ADC 
and the DSP on chip enabling high quality detection com- 
bined with incomparable portability (El Gamal et al, 1999). 
Intelligent use of CMOS circuitry can be used to compensate 
for non-idealities engendered by the poorer quality (relative 
to CCDs) of added noise from the readout chain. Inclusion 
of transparent on-chip background subtraction and signal 
averaging circuitry are two examples that can help improve 
CMOS detection quality beyond that of CCD's. Background 
subtraction removes the deterministic thermally generated 
portion of the photodiode signal, leaving only its shot noise 
component that has a variance equal to the square root of its 
mean. Signal averaging lowers the independent noise com- 
ponents by yfn, where n is the number of independent sam- 
ples. Indeed, using CMOS as a biological platform has the 
potential of making "in the field" biological testing a reality. 

A conceptual schematic of the proposed integrated sys- 
tem in which the chemistry is directly coupled to the 
CMOS-based detection chip is shown in Fig. 5a. A block 
diagram of the proposed detection chip is shown in Fig. 5b. 
It consists of a 2D array of pixels, each containing a pho- 
todetector and transistors for readout. The collected charge 
from each pixel is read out and converted in parallel to a 
digital format via an array of per-pixel analog-to-digital con- 
verters (ADCs). The system includes memory for on-chip 
storage of multiple frames. Moreover, a dedicated digital 
signal processor (DSP) is integrated on the chip to perform 
any needed computations, such as for background subtrac- 
tion and read noise reduction. A detailed description of a 
prototype is provided in (Eltoukhy et al., 2004). 

In Sections 3.1, 3.2 and 3.3 we describe the modeling 
and algorithms involved in the optics, image sensor and 
post-processing components of the detection chain. 

3. 1. Optics model 

We first compare the optical efficiency of direct coupling 
to a conventional camera-based imaging system. Assuming 
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Fig. 5. (a) Conceptual schematic of the proposed integrated platform demonstrating the tight coupling between the sensor array and the biological assay 
reactions and (b) block diagram of the detection chip showing its main components. 
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that the object to be imaged is a point source located a 
distance U away from the lens (or the imaging surface in 
the case of direct coupling), the numerical aperture NA of 
the object side can be approximated using geometric optics 
by sin(<9 max ), where 0 max — tan~ x (D/2U) where D is the 
diameter of the lens aperture. Since optical efficiency (i.e. 
fraction of light collected by the lens) is proportional to 
(NA) 2 , the efficiency of a camera-based system (see Fig. 6a) 
for typical parameter values of D = 0.5 f and U = 2f, 
where/ the focal length is a mere 1.6%. Optical efficiency 
can be greatly improved by directly coupling the chemistry 
to the detection chip surface as proposed in our integrated 
system. For example, taking D = 1 00 jxm (photodetector 
size) and U = 15 firn, the optical efficiency becomes 91%. 

To accurately compute the optical efficiency for the 
direct coupling scenario, we use an optical simulator 
(LightTools™) and the specific luminescent assay and sen- 
sor parameters of the proposed setup to compute the point 
spread function, PSF (i.e. the function describing the distri- 
bution of light at the imaging plane due to a point source at 
the object plane). The PSF relates to the degree of blurring 
encountered in an optical system and can be used to calculate 
the resultant light intensity at the image plane due to the finite 
separation distance between the assay and the sensor. This is 
performed by convolving the derived PSF with the intensity 
distribution at the object plane to obtain the corresponding 
distribution at the imaging plane. Fig. 6b depicts a simpli- 
fied model of the imaging setup and the simulated intensity 
distribution at the sensor plane is plotted in Figs. 6c and d. 



3.2. Image sensor model 

An image sensor comprises an array of pixels each hav- 
ing a photodetector and devices for readout. The sensor 
is typically operated in direct integration as illustrated in 
Fig. 7. Photons incident on the photodetector are converted 
into photocurrent. The photocurrent is directly integrated 
over the photodiode capacitance Cd into charge Q. The 
amount of charge that can be collected is limited by the 
well capacity 2max- At the end of integration time T\ nU the 
charge is read out as a voltage signal V 0i which is related 
to Q by the conversion gain g, i.e. Vb = gQ- Fig. 8 sum- 
marizes the overall image sensor model, including dark 
current, added shot and read noise, and the conversion 
gain. 

The efficiency of converting incident photons to photocur- 
rent / p h(0 is governed by the spectral response r)(k) of the 
detector, where k is the emission wavelength, and is given 
by 

*ph(0 = r](X)I inc[dcnt (t). (20) 

Several sources contribute to noise during the collec- 
tion of the photogenerated signal including dark current 
z'dc and photocurrent shot noise, reset noise, and readout 
noise. The total added average noise power can be expressed 
as 

o2 = £(Q + Gdc) + eLd. < 21 ) 
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Fig. 8. Image sensor model, including added dark current, shot and read noises. 



where the first term is the average power of the integrated 
shot noise of both the signal and the dark current, Q 2 Qad is 
the average read noise power, and q = 1.602 x 10~ 19 C. 
The signal-to-noise ratio, SNR, is defined as the ratio of the 
photogenerated signal power to the noise power and is given 
by 

Q 2 

SNR =^=-. (22) 

°n 

3.3. Image post-processing 



the optics and image sensor. For instance, the PSF can be 
used to correct for assay cross-talk using conventional equal- 
ization techniques, such as Weiner filtering. This enables 
higher throughput for the detection system, since tighter as- 
say pitches can be tolerated. Fig. 9 shows an example of this 
filtering technique on an image with excessive cross-talk. 
Additionally, subtraction of both the chemical and detector 
backgrounds as well as digital accumulation incorporating 
the reaction kinetics can be performed to enhance SNR and, 
thus, detection reliability. 



Once the image data is collected, post-processing algo- 
rithms can be applied to correct for non-idealities in both 



= !■ . «■ <. v, ■ s | 

w 




(b) 



Fig. 9. Microarray image (a) before post-processing and (b) after applying 
cross-talk reduction algorithm. 



4. Simulation and experimental results 

By combining the models and algorithms discussed 
in the previous sections for the reaction kinetics, opti- 
cal path, image sensor and post-processing, the minimum 
detectable analyte concentration can be estimated for a 
given set of detection system parameters. A block dia- 
gram of the luminescence detection system simulator that 
implements these models is shown in Fig. 10. The sim- 
ulator, including the reaction kinetics and image sensor 
models, are coded in Matlab, while the optical efficiency 
and PSF for a given system geometry are calculated us- 
ing LightTools™ and passed as arguments to the Matlab 
script. 

To validate the above simulation model, we built an 
isolated imaging chamber for luminescence detection. The 
system employs an ultrasonic sprayer attached to a robotic 
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kt : Turnover rate of ATP 
ka : Turnover rate of Luciferase 
a : Luciferase Efficiency 
V : DNA Volume 



System Geometry 
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Tj(X) : Spectral response 
Idc : Dark Current 
Tint: Integration Time 
G r : Read Noise 



Fig. 10. Block diagram of luminescence detection system simulator. 



Background subtraction 
Cross talk reduction 
SNR Enhancement 



arm for nucleotide and enzyme delivery to the target assays 
above the imaging apparatus as shown in Fig. 1 1 . In order to 
simulate the types of conditions that would be encountered 
in an integrated CMOS detection system, a scientific-grade 
CCD (Hamamatsu HC230) is incorporated with the ultra- 
sonic sprayer. The system allows direct placement of the 
assay slide onto the CCD. The CCD output is connected to 
a dedicated PC for simultaneously controlling the robotic 
arm of the sprayer, capturing images, and processing the 
assay data to achieve accurate detection. We implemented 
a simple but efficient algorithm for signal detection that 
can run in parallel with the data acquisition software. The 
algorithm automatically performs background subtraction 
of noise and utilizes line binning to increase SNR. An 
eight-tap low pass filter is utilized to smooth the resultant 
signal and remove any high frequency noise components. 
Also integrated with the script is an equalization routine for 
reduction of cross-talk using the modeled system PSF (see 
Fig. 12). This experimental setup allows for accurate mea- 
surement of the quantum efficiency of the reaction as well 
as of the optical and detection paths. The overall system 
dimensions are 50 cm x 50 cm x 60 cm due to the pres- 
ence of the robotic arm. In contrast, the envisioned CMOS 
system under development will measure 10 cm x 10 cm 
x 2 cm. 



We have used the experimental system to measure the 
light intensity produced during the incorporation of a sin- 
gle base in a Pyrosequencing reaction. Using lOOfmol of 
DNA per well, which is 10-fold lower than the amount 
currently used in commercial Pyrosequencing systems, we 
measured a peak signal-to-noise ratio of 750. We have also 
used the collected experimental data to fine-tune the param- 
eters of the chemical kinetics, optics, and detector models. 
This "tuning" involved incorporating the specific associa- 
tion and disassociation rate constants, optical path loss, sen- 
sor noise characteristics, and overall system delay resulting 
from the finite mixing time of the sprayer. Fig. 13 com- 
pares the simulated and experimental results during incorpo- 
ration of one nucleotide in a Pyrosequencing reaction using 
lOOfmol of DNA. As is readily seen, the simulation results 
are fairly well corroborated by the experimental data. Fur- 
thermore, using typical CMOS process parameter values at 
room temperature, our model predicts that 0.5 photon/s/fim 2 
in conjunction with 7j nt = 10s and 100 \xm x 100 [xm diode 
sizes would be needed to achieve a signal-to-noise ratio of 
10. This implies that an integrated CMOS-based lumines- 
cence detection system would be capable of performing py- 
rosequencing with as little as 1 fmol of DNA, an amount 
three orders of magnitude lower than current commercial 
machines. 




Fig. 1 1 . Image of the prototype system including a robot-arm, an ultrasonic sprayer and a cooled CCD sensor. 
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Fig. 12. (a) Plot of measured output signal voltage from each column of the CCD array for two adjacent 0.1 pmol wells and (b) same output after 
applying signal processing algorithms to reduce the effect of the system PSF. 
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Fig. 13. Experimental data vs. combined chemistry and sensor simulation 
of detected signal vs. time of nucleotide incorporation in a Pyro sequencing 
reaction. 



5. Conclusion 

We described a simulation model for CCD and CMOS- 
based luminescence detection platforms. The model quanti- 
fies the photon flux generated by luminescence probes using 
ATP-based and luciferase label-based assay kinetics. The 
photon flux coupled with the system geometry is then used 
to calculate the incident photon flux on the imaging plane. 
Subsequently, the output image is computed using a de- 



tailed image sensor model. We constructed a prototype sys- 
tem to experimentally verify the developed models. Using 
this experimental setup we are able to obtain accurate mea- 
surements of the quantum efficiency, temporal kinetics and 
spatial distribution of the reaction, which are used to cal- 
culate the optimal assay sizes and throughput limits for the 
CMOS-based system. 
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